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ABSTRACT 

Recently Sanwal et al. (2002) reported the first clear detection of absorption features in an isolated 
neutron star, 1E1207. 4-5209. Remarkably their spectral modeling demonstrates that the atmosphere 
cannot be Hydrogen. They speculated that the neutron star atmosphere is indicative of ionized Helium 
in an ultra-strong 1.5 x 10^** C) magnetic field. We have applied our recently developed atomic 
model (Mori &: Hailey 2002a) for strongly-magnetized neutron star atmospheres to this problem. We 
find that this model, along with some simple atomic physics arguments, severely constrains the possible 
composition of the atmosphere. In particular we find that the absorption features are naturally associated 
with He-like Oxygen or Neon in a magnetic field of ^ 10^^ C, comparable to the magnetic field derived 
from the spin parameters of the neutron star. This interpretation is consistent with the relative line 
strengths and widths and is robust. Our model predicts possible substructure in the spectral features, 
which has now been reported by XMM-Newton (Mereghetti et al. 2002). However we show the Mereghetti 
et al. claim that the atmosphere is Iron or some comparable high-Z element at ^ 10^^ C is easily ruled 
out by the Chandra and XMM-Newton data. 

Subject headings: stars: neutron ~ X-rays: stars - stars: atmospheres - individual: 1E1207. 4-5209 



1. INTRODUCTION 

A major goal of neutron star (NS) research has remained 
unrealized despite 30 years of effort - to determine the fun- 
damental properties of the superdense matter in a NS inte- 
rior, in particular its equation of state (EOS) and compo- 
sition. There are many approaches which can be employed 
in this effort (van Kerkwijk (2001) and references therein), 
but one of the most intensively studied is to exploit the 
thermal radiation from isolated NS. This radiation can be 
used to deduce information about the EOS from neutron 
star cooling theory (Tsuruta et al. 2002; Yakovlev et al. 
1999). Extracting information about the NS interior is 
not straightforward, however, since the observed spectrum 
does not represent the NS surface emission, but is mod- 
ified by radiative transfer effects in the NS atmosphere. 
The problem of unfolding the observed spectrum and un- 
derstanding NS interiors thus depends on deducing the 
composition of the atmosphere. 

NS atmospheres are interesting in their own right. 
We know little about their nature, other than that 
some NS probably have Hydrogen atmospheres (Pavlov 
et al. 2002a), since they fit well to the sophisticated H- 
atmosphere models that have been developed (Pavlov et al. 
1995; Zavlin & Pavlov 2002). The observations have been 
silent on the question of whether non-Hydrogen atmo- 
spheres exist. 

The NS 1E1207.4-5209 was discovered by Helfand & 
Becker (1984) with the Einstein Observatory and is associ- 
ated with the SNR PKS 1209-51. Subsequent observations 
established it as a prototypical radio silent NS associated 
with a SNR (Bignami et al. 1992; Caraveo et al. 1996). 
ROSAT and ASCA observations were fit with a black-body 
spectrum (Mereghetti et al. 1996; Vasisht et al. 1997) and 
a Hydrogen atmosphere model (Zavlin et al. 1998). No 
X-ray pulsations were detected. Chandra detected X-ray 



pulsations with P = 0.424 s (Zavfin et al. 2000). While 
this is perhaps not surprising in light of previous work, 
the small period derivative estimated (Pavlov et al. 2002b) 
certainly was unexpected. The inferred surface B-field is 
- (2 - 4) X 1012 G. 

Even more remarkable, Sanwal et al. (2002) (hereafter 
SZPT) have discovered features at ~ 0.7 keV and ^ 1.4 
keV are required to obtain acceptable fits to the spectrum 
(which was fit with a 2.6 x 10^ K underlying black-body 
continuum) . A feature of marginal significance and unclear 
origin was also noted at 2 keV. A subsequent observa- 
tion by XMM-Newton (Mereghetti et al. 2002) has pro- 
vided some hint of substructure in the features, and has 
marginally detected the ^ 2 keV feature seen by Chandra. 

Three different phenomenological fits indicative of ei- 
ther absorption lines or edges were used by SZPT to fit 
the two strong features. Of enormous significance, the 
spectrum was shown to be inconsistent with Hydrogen 
atmosphere models. The intense effort which has gone 
into the Hydrogen atmosphere models and their success 
in explaining other NS observations (Pavlov et al. 2002a) 
lends credence to the inevitable conclusion that the at- 
mosphere of 1E1207. 4-5209 is something other than Hy- 
drogen. Unfortunately, as pointed out elsewhere (Zavlin 
& Pavlov 2002), work on non-Hydrogen atmospheres at 
high B-fields is much less developed. SZPT argued on 
various grounds that the features could not arise via cy- 
clotron lines. Instead they tentatively suggested emission 
from a once ionized-Helium atmosphere with a B-field of 
1.5 x 10^'* G. This B-field is inconsistent with that derived 
from the spin parameters, but SZPT argue this could be 
due to an off-centered B-field or glitches affecting the P 
measurement. Others have argued for a cyclotron line so- 
lution at lower B-field (Xu et al. 2002). Mereghetti et al. 
(2002) claimed an atmosphere of Iron or other high-Z el- 
ements at a B-field of ~ 10^^ G, although they did not 
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actually fit their XMM-Newton data to a model. 

In this paper we offer alternate interpretations of these 
spectral features. All our interpretations involve atomic 
transitions, mainly in He-like ions of mid-Z elements at 
B ~ 10^^ G, consistent with the B-ficld derived from the 
NS spin properties. The most likely of these interpreta- 
tions is that the NS atmosphere contains Oxygen or Neon; 
most noteworthy is that all our models, whether consid- 
ering just the two strong features or including the third 
weak feature, demand mid-Z elements for an acceptable 
solution. Our model, combined with the XMM-Newton 
and Chandra data, easily rule out the Iron and high-Z so- 
lutions of Mereghetti et al. (2002). 

Some comments are in order on our approach. The 
atomic spectroscopy data used in this analysis is based 
on a novel approach for obtaining fast and accurate so- 
lutions to the Schrodinger equation for B-fields in the 
Landau regime (appropriate for all cases considered here). 
This approach, multiconfigurational, perturbative, hybrid, 
Hartree, Hartree-Fock theory, allows rapid computation 
of transition energies and oscillator strengths for arbi- 
trary atom, ion, excitation state and B-field (Mori & Hai- 
ley (2002a), hereafter MH02a). This permits a complete 
search of all possible spectroscopic transitions consistent 
with the given line or edge energies. 

While it may appear that this approach produces 
an uninterestingly large number of potential solutions, 
we demonstrate in a companion paper (Mori & Hailcy 
(2002b), hereafter MH02b) that this is not the case. We 
show that the presence of two or more line or edge fea- 
tures provides a remarkable robustness to a host of poorly- 
understood atomic physics effects and unambiguously re- 
stricts the atmosphere composition to mid-Z elements. 

2. Chandra data reduction 

We only briefly mention our data reduction here, as our 
approach closely follows that of SZPT. Indeed, we empha- 
size that none of our conclusions here would be modified 
if we simply used the spectral line parameters derived by 
SZPT. Our analysis is described in more detail in MH02b. 
The results presented here are for phase-integrated spectra 
only. Subsequent work will consider phase-resolved data, 
which can provide more information on the system geom- 
etry. We fit the spectrum with two models: a blackbody 
with two absorption edges (edge model) and a blackbody 
with two absorption lines (line model). As noted by SZPT, 
these spectral features arc required to get acceptable fits. 
The fit parameters for our models are given in Table 1. 
We show the raw Chandra spectrum in figure 1 overlaid 
with a blackbody to clearly indicate the presence of the 
absorption features. Figures 2 show the results of our best 
fit line and edge models. 

3. GENERAL SPECTROSCOPIC CONSIDERATIONS 

We briefly review the quantum numbers of bound elec- 
trons in B-fields in the Landau regime. States are labeled 
by (m, f) where m is the magnetic quantum number and 
v is the longitudinal quantum number. The v = Q states 
have larger binding energy than the > states. Consis- 
tent with the terminology of MH02a, we call the former 
tightly-bound states and the latter loosely-bound states. 
Generally the ground state of an ion is a tightly-bound 



state. In the X-ray band transitions are mostly likely to 
be either photo-absorption (tightly-bound state to contin- 
uum) or absorption lines from tightly-bound to loosely- 
bound state (tight-loose transition). Tight-tight transi- 
tions appear in the optical band. Details of the transition 
selection rules are discussed elsewhere (MH02a,MH02b). 

4. specific atmosphere solutions (two spectral 

features) 

We first considered the two strongest features at ~ 0.7 
keV and ~ 1.4 keV and assume both features are due to 
atomic transitions in a single element and ionization state 
(Case A). This effectively defines a solution for the domi- 
nant absorbing element in terms of He- like ions. Li- like 
ions and those in lower ionization states are not dom- 
inant absorbers under the expected atmospheric condi- 
tions (MH02b). H-like ions cannot be the dominant ab- 
sorbers because this is grossly incommensurate with the 
features' relative widths and strengths even under condi- 
tions far removed from local thermodynamic equilibrium 
(MH02b). H-like ions can, however, be present since the 
relevant spectral features will be blended with the He-like 
features. In the He-like ions, the electrons are making 
tight-loose transitions out of the (0,0) (1,0) ground state 
to the (0, 1)(1, 1) excited states respectively, or to the con- 
tinuum. Neither tightly-bound states at higher magnetic 
quantum number nor non-adjacent tightly-bound states 
are acceptable. In the former case the feature energy ra- 
tios cannot be produced and in the latter case unobserved 
features in the Chandra band would be produced. 

We searched for any combination of element, redshift 
and B-field which could give the observed positions of 
the spectral features. Our B-field range was ~ 10^^ G 
to ~ 2 X 10^'' G. In the case of edge solutions we also per- 
mitted (small) pressure shifts (details in MH02b). Line 
solutions are insensitive to pressure shifts. The properties 
of the solutions including derived values of the B-field and 
redshift are shown in table 2. The first two solutions in- 
volve atomic line transitions in He-like Oxygen and Neon 
and the third solution He-like Neon absorption edges. A 
similar Oxygen solution was rejected because its redshift 
was too low, violating causality and stability conditions 
for the NS interior (Lindblom 1984; Haensel et al. 1999). 
For the Neon edge solution there is an electron cyclotron 
line in the Chandra/ XMM-Newton energy band but it was 
made to overlap the higher energy spectral feature by in- 
troducing a small pressure shift (^ 50 — 90 eV). Such 
pressure shifts slightly affect the energy ratio of the two 
edge features and serve to further constrain other elements 
and conditions of the atmosphere (MH02b). The line and 
edge widths are all consistent with expectations for Stark 
broadening and non-constant B-field broadening respec- 
tively (MH02b). 

Our case B assumes that one spectral feature is an 
atomic transition and the second is a cyclotron line. A 
1.4 keV or 0.7 keV electron cyclotron line defines a B-field 
of 1.2 xlOii(l-F^) G or 0.6x10^^1-1-0) G respectively. In 
the former case solutions can only be found for H-like ions 
transitioning out of the (0, 0) ground state to the (0, 1) ex- 
cited state or the continuum. These are all mid-Z solutions 
albeit with lower B-field than those considered above. In 
the case of a 0.7 keV electron cyclotron line, the B-field is 
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low enough that our perturbation method breaks down for 
higher Landau levels and is therefore loss accurate; we es- 
timate that putative solutions would involve the elements 
Beryllium or Boron and are therefore implausible. If one of 
the features is a nuclear cyclotron line then the other fea- 
ture must be a transition from the ground state of either 
Helium or Hydrogen. However there is no combination 
of B and z that can reproduce the line positions. Thus 
nuclear cyclotron lines are ruled out. 

Although Case B is formally allowed wc consider it 
highly unlikely. There are difficulties in attempting to 
reconcile the depth and width of the spectral features 
(MH02b) . Nevertheless it is interesting that even this sce- 
nario including cyclotron lines only admits solutions in- 
volving mid-Z elements. 

5. specific atmosphere models (three spectral 
features) 

Table 2 also shows proposed solutions if the feature at 
~ 2 keV is confirmed (Case C solutions). Again it is of 
importance to note that only mid-Z solutions are accept- 
able. We had to relax one of our previous assumptions 
by assuming additional significant absorption from H-like 
ions. An admixture of H-likc ions in a given atmosphere 
is qiiite plausible for reasonable conditions of temperature 
and density (MH02b). 

The first case considered is identical to the He-like Oxy- 
gen/Neon case presented above but with the addition of 
a H-likc Oxygen/Neon absorption edge due to the (0,0) 
ground state to continuum transition. In this model the 
other H-like lines overlap the He-like lines so no additional 
features arc produced. 

The second solution invokes photoabsorption edges in 
He-like Neon for the two lower energy features and the 
higher energy feature is due to an electron cyclotron line. 

6. discussion 

Our most important conclusion is that all viable solu- 
tions for B-fields comparable to that inferred from the NS 
spin parameters require mid-Z atmospheres in both the; 
two and three spectral feature case. Our solutions which 
most closely correspond to the inferred B-field demand 
He-like Oxygen and Neon. 

In contrast to this paper, SZPT proposed ultrahigh B- 
field solutions. Their Helium atmosphere appears prob- 
lematic, however. In order to get the right feature po- 
sitions SZPT's H-like Helium features must of necessity 
arise in transitions from electrons in the (0,0) or (1,0) 
quantum states to (0, 1) and (1, 1) states. The (0, 1) state 
is heavily depopulated at the types of ion densities one ex- 
pects and the (1,1) state is actually autoionized due to the 
large nuclear cyclotron energy. In fact the motional Stark 
shift and pressure effects in Helium will likely destroy all 
the loosely-bound states. Consequently one is forced to 
assume the features are absorption edges, and this leads 
to severe difficulty explaining the relative edge strengths 
(MH02b). 

Mcrcghetti et al. (2002) did not fit their XMM-Newton 
data, but proposed nevertheless that the spectral features 
could be attributed to Iron or some other high-Z element 



at a field ~ 10^^ G. Our current work rules out their pro- 
posal. There is no combination of B-field, redshift, element 
or ionization states, even with unphysical, arbitrary pres- 
sure shifts, that is consistent with the XMM-Newton and 
Chandra data. 

The interpretation we present here makes an important 
prediction: the observed features may show substantial 
substructure when observed with higher spectral resolu- 
tion. Transition energies from the tightly bound ground 
state to loosely bound final excited states differing in lon- 
gitudinal quantum number are rather close to each other, 
and may be blended. Details of predicted energies and 
transition assignments are presented in MH02b. Whether 
detailed atmospheric models would lead to a washing out 
of this substructure through broadening is not clear, but 
its detection would certainly confirm that the features are 
atomic transitions since we would not expect similar sub- 
structure in cyclotron lines. In fact, Mcrcghetti et al. 
(2002) claim a marginal detection of such substructure. 

Concerning the B-field derived from the free parameter 
fits to the data, we note that our Neon and Oxygen line so- 
lutions are all within a factor of two or three of the B-field 
derived from spin parameters. The discrepancy may lie in 
the NS B-field geometry or line blending. Line blending 
will affect our B-field (and redshift) estimate, but will not 
affect our element identifications. 

The analysis here, along with that of SZPT, begs the 
question of why there is anything in the atmosphere other 
than Hydrogen. Oxygen or Neon are expected constituents 
near the NS mass cut in many models of supernova explo- 
sions (Arnett 1996). Indeed due to turbulent mixing in 
such explosions Oxygen can be found throughout the core 
(Herant & Woosley 1994). But given that a very thin 
layer of Hydrogen accreted from the interstellar medium 
is enough to produce an optically thick Hydrogen atmo- 
sphere, wc are left to ponder mechanisms by which such 
accretion can be inhibited, or at least permit simultane- 
ous emission from a mid-Z element. One possibility is that 
the timescale for an optically thick layer to accrete onto 
the surface is shorter than that for gravitational stratifica- 
tion in the atmosphere (van Kcrkwijk 2001). In that case, 
even trace amounts of mid-Z elements could lead to strong 
emission (Pavlov et al. 2002a). Optical observations of the 
associated SNR indicate the presence of Oxygen and it has 
been argued that the unusually high Galactic latitude of 
this NS implies the Oxygen originated in the progenitor 
star (Ruiz 1983). Thus a rich source of Oxygen is available. 
Alternately the mid-Z atmosphere may be associated with 
the NS mass cut, and a propeller effect may be severely 
inhibiting accretion of Hydrogen. We note that several of 
our proposed solutions allow a Hydrogen-rich environment 
(MH02b). 

For more progress to be made on 1E1207. 4-5209 it will 

be necessary to develop complete models of NS atmosphere 
for mid-Z elements incorporating the more sophisticated 
atomic data bases we have utilized here. This source has 
much to tell us in the years to come, especially when higher 
resolution spectroscopic observations become available. 

The authors thank David Helfand for helpful discussions 
and a careful reading of the manuscript. 
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Table 1 

Results from spectral fitting by the three models. 



Edge Gaussian Lorentzian 



kT°° [cV] 


269 ±2 


254 ±2 


254 ±2 


Nh [102"cni-2] 


4.5 ±0.3 


4.6 ±0.3 


5.1 ± 0.4 


[km] ^ 


1.68 ±0.04 


1.80 ±0.05 


1.91 ±0.07 


El [eV] 


626 ±5 


737 ±6 


737 ±6 




0.60 ± 0.05 


0.46 ± 0.03 


0.57 ±0.03 


EWi [cV] 


338 ± 28 


119 ± 10 


327 ±27 


E2 [cV] 


1298 ± 10 


1424 ± 10 


1430 ± 10 




0.54 ±0.05 


0.47 ± 0.04 


0.55 ± 0.05 


EW2 [oV] 


637 ±53 


115 ±10 


199 ± 17 


line (edge) ratio 


2.07 ± 0.03 


1.93 ± 0.03 


1.94 ±0.03 


^/ 


1.07 


1.01 


1.10 



'^We assumed the distance ~ 2.1 kpc for apparent radius. 
''Optical depth of edges and lines. 
■^142 degrees of freedom. 



Table 2 

Candidates for elements, magnetic field strengths and gravitational redshifts. 



Case 


Ion 


Features 


B12 


Redshift 


Pressure shift? 


A 


VII 


two lines 


0.55-0.75 


0.06-0.21 


No 


A 


Ne IX 


two lines 


0.80-1.10 


0.62-0.86 


No 


A 


Ne IX 


two edges 


0.18-0.20 


0.56-0.57 


Yes 


B 


C VI 


one edge & cyclotron line 


0.21-0.22 


0.70-0.73 


No 


B 


N VII 


one line & cyclotron line 


0.13 0.14 


0.10-0.13 


No 


B 


VIII 


one line & cyclotron line 


0.18-0.19 


0.45-0.49 


No 


C 


VII, VIII 


two lines & one edge 


0.55-0.75 


0.06-0.21 


No 


C 


Ne IX, Ne X 


two lines & one edge 


0.80-1.10 


0.62-0.86 


No 


C 


Ne IX 


two edges & cyclotron line 


0.26-0.28 


0.59-0.60 


Yes 



